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Abstract 

The holographic dual of self gravitating Julia-Zee dyon is discussed. It is shown that the 
dual field theory is generally a field theory with a vortex condensate. The vacuum expectation 
values of the dual operators, as functions of the sources in the field theory, are studied in a 
class of bulk solutions. In these solutions the sign of the vacuum expectation values of the dual 
operators change, by changing the sources in the model. 



1 Introduction 



It is well known that the AdS/CFT correspondence [l]-[3], is a useful tool in understanding at least 
the qualitative features of strongly coupled field theories. During past years many holographic 
models are studied. In some cases the bulk theories are constructed so as to simulate certain 
features of a known quantum field theory, e.g., AdS/QCD and AdS/CMT models; and in some 
cases a field theory is proposed as the dual theory of a given bulk theory. In both approaches, 
by the well known recipe of AdS/CFT correspondence, the key considerations are symmetries and 
the field contents of the theories. In some cases the comparable quantities and behaviors are 
illuminating for finding the dual theories. 

In this work we study one of the gravitational systems recently studied, the self gravitating Julia- 
Zee dyon [4j|6]. A.R.Lugo, E.F.Moreno and F.A.Schaponsik, have found a numerical solutions for 
the Julia-Zee black dyon, in the BPS limit, with [s] and without [4] back reaction on the metric, 
and they found that the vacuum expectation values vanish by raising the temperature of the bulk 
blackhole. It can be considered as a phase transition from a dyonic blackhole. In [g], it is shown 
that a simple, vanishing scalar solution for the dyon, becomes dynamically unstable by lowering 
the temperature of the bulk blackhole to a scalar hairy blackhole, which means that the dual field 
theory changes its vacuum, so by lowering the temperature of the dual boundary system, a phase 
transition occurs to a scalar condensed phase. Multimonopole solutions with AdS background are 
studied with axial (7|[8] and crystalline [9, 10 symmetry. In [9], D.Tong and S.Bolognesi have 



studied the multi monopole solution at the BPS limit, and found that at low temperatures, the 
multi monopole solutions are favored in comparison to blackhole single monopole solutions. As the 
multi monopole solutions in AdS background can have crystalline symmetry, they have proposed 
that these configurations could correspond to strange metals at low temperatures. In this article 
we study the holographic vacuum expectation values for the dual operators, by the well known 
recipe of the Gauge/Gravity duality; and also calculate a numerical solution to the equations of 
motion, to determine the behavior of the v.e.vs of the dual operators, as functions of the sources. 

The paper is organized as follows. In the next section we introduce the bulk action, equations 
of motion and the form of the Julia-Zee dyon. In section (3) we calculate the vacuum expectation 
values of the dual operators, generally from Einstein- Yang-Mills-Higgs action and as a special case, 
for Julia-Zee dyon. In section (4) we find a set of numerical solutions for the dyon, and generically 
study their thermodynamics and v.e.v of the dual operators. In the last section we end with 
conclusion. 

2 The model 

The action we consider in this article is 

s =Jv^gd'x [^[R+ ^) - If^^f'^^'^ - iD.cp^D^r - ^{rrf + ^irr)], (2.1) 

The covariant derivative and the field strength are: 

D.cp- = d.^P"" + ee^'^A'^^' ; F^, = d^Al - d.A^ + ee^'^A^Al. (2.2) 

As is clear, the mass term for the scalar field is chosen to be m? = —2/ LP'. Since the su(2) gauge 
symmetry of the bulk action corresponds to a global su(2) symmetry in dual field theory, so the 
scalar field can be considered as the dual of a two fermion flavored composite operator (similar 
to pions or a flavored Cooper pair). In this sense, since the mass dimension of a fermion in 2+1 
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dimensions is one, the operator dual to the (j), in the field theory, will have mass dimension two. 
On the other hand the source that couples to such an operator has mass dimension on^ In order 
to have these special mass dimensions, the asymptotic form of the bulk scalar field should be: 



~r 9 ~r ••• • 



(2.3) 



So vr? , has to be —2/L'^. 

The equations of motion are: 



1 6 

-{R+j^)g^u = 8TrGT^, 



D„ 



5V 



(2.4) 



-g e 



abc lb ( T^v i\c 



The Julia-Zee ansatz has the form: 

- H(r) 
er 



A = a.A = ^ a' dt + ^ — {- dO + sin{e) dip), 



(2.5) 



er 



where ,cr^ , o'-p are defined as: 



a = r.a , a 



9.a , a"^ = (2.6) 

and f , 6 , ip are the usual unit vectors of spherical coordinates We consider the metric in the 
most general form of a spherical symmetric metric in global coordinates: 



ds' 



There are two simple solutions for the full equations of motion [6,11 



H = Cr , K = 0, J = fir-p, = e'^ = 1 - — + % + 



2 ^2 



L2 



H = , K = , J = ^ir-p , e 



X Y , 2M q 



2 j,2 



L2 



where C = is the minimum of the potential, and q is related to p through: 

2 _47tGn{1 + P^) 



(2.7) 



(2.8) 



(2.9) 



and p is related to the chemical potential by 



P = prn- 



^It is similar to a Dirac or Majorana mass, 
■^we name this frame as spherical frame of su(2). 



(2.10) 
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L satisfies: 

It is very easy to see that for the above solution there is a relation between the temperature of the 
blackhole and chemical potential as: 

T + a^'^ = l3 (2.12) 

where a and j3 are two constants related to the parameters of the solution {X,rH,GN,e,L). This 
relation is in fact the first law of thermodynamics: 

=T + u - ; poc Li (2.13) 
s s 

3 Vacuum of dual theory 

To understand the vacuum of the dual field theory, we should calculate the vacuum expectation 
values of the operators in the field theory. In the theory under consideration we have a scalar field, 

which is dual to a scalar operator, O", and a gauge field A^, which corresponds to a vector 
operator, V^, in the dual field theory. 

Before we proceed to calculate the one point functions of the operators, we recall the linear 
method for calculating them in AdS/CFT. 

Suppose that we have a bulk theory which is defined by an action S(^). In order to calculate 
the 1-point functions we need to vary the on-shell action with respect to the fields. Expanding to 
the second order we have: 

Son-sheiiicl^B + 5(t>) = S^^\cI)b) + 5«(0b,(5</.) + S^''\4>B,5<t>'') + (3.1) 

the first term is the free energy of the dual theory. The second and third terms can be written in 
the forms: 



= j (5(/).(b.e.m) + S^^\(l)B,m 
5(2) = j 5<^.(f.e.m)+ S^^\<t>B.5ct>^) 



(3.2) 



the terms b.e.m and f.e.m refer to the expressions for the background equation of motion and that 
of the linearized equation of motion for the fiuctuation, respectively; which are zero. Now suppose 
that the asymptotic expansion of the fields are: 



4>B = 4>sr^+ + ■■■ + 4>cr^- + 
d(j> = S^pr^+ + ... + S4>rr^- + 

where the subscripts s,c,p,r, refer to source, condensate, probe and response, respectively. 
In general the relevant parts of the Sn's are of the forms: 

sj;^ = / ... + CB 0c + - 



(3.3) 



= J ... + Cr 54>p 54>r + 



(3.4) 



where cs are constants. Then the one point function for the operator dual to (p will be: 

(O) = =CB^c + Cr S4>r = {O) B + {0)r, (3.5) 
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where subscripts B and r refer to the background contribution to the v.e.v of O, and the response 
of the system to the probe, respectively. 

Doing the above calculation for the Einstein- Yang-Mills-Higgs action leads to: 

S{Ab + 6 A, (^B + S(l)) =S^''\Ab, (t>B)+ 



sf{AB, (PB, S^^) + S^^liAs, ^B, SA, 6(t>) + S^^\Ab, ^b, SA"") + 



(3.6) 



=S^'\AbAb)+ 
s'iliAB, ^B, H) + 4a (^B. <t>B, 5A)+ 
S^^\{Ab, <Pb, 50') + S^^\{Ab. ^b. 5A^) + 

In the last line we have used the equations of motion for background and linear equations for 
fluctuations. For Einstein- Yang-Mills-Higgs action, the Noether actions are: 

Sn\ = - [ d^.D^'ch ds^ = - 6$.D'cbB d^x 



Sn\ = -\j ^9 5^.D^^5$ ds^ = -lj^g 6^.D^6$d^x 
= -JV^ dA^.F^" ds^ = - 6A,.F^'' d^x 



(3.7) 



S = -11 5A,.5F>^^ ds^ = Al-r-.>^A,T?ru ,3. 



where ds^ is the element normal to the boundary. 
So roughljj^we have: 

{V^)b CC (FB)rM , {d)B OC OJb- (3.8) 

The v.e.v of the charge density operator (Vq = p), is proportional to the radial electric field; 
the v.e.v of spatial components of the vector operator (Vi), are proportional to the tangent to the 
boundary components of the magnetic field, and the v.e.v of the scalar operator is proportional to 
the radial derivative of the scalar field. So for Enstein- Yang-Mills-Higgs action, the radial magnetic 
field, which is the characteristic of magnetic monopoles, doesn't contribute to v.e.vs as a condensate; 
but, its boundary value plays the role of a source for the dual field theory. If we want the radial 
magnetic field to contribute to the v.e.vs, as a condensate, we have to add other terms to the action. 
For example adding a 0-term leads to: 

(p) OcEr + d Br 

: ^ ^ (3.9) 

{Vi) OC Cij {Bj + e Ej) 

3.1 Julia-Zee dyon and Vortex condensation 

Applied to the Julia-Zee dyon, it is easy to see that: 

{Vt) « {E)r OC (^)' f 

^ ^ r 

(Vi) OC e-ij{B)j OC K' eijdjf 

■^By "oc", we mean the v.e.vs are proportional to one of the expansion coefficients which couples to the probe, as 
mentioned previously. 
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(a) (b) 
Figure 1: plot of directed curves for e+ (a) and e_ (b) basis. 

(0) oc 9^0oc (— )' f. (3.10) 
r 

Thus for both solutions in (2.7), these quantities are zero, except for the charge density. In the 
next section we will introduce a numerical set of solutions for which the vacuum expectation value 
of the operators are non zero, and we will study them in the phase space; but, before we proceed, 
it is interesting to look at the profile of the vector condensate. 

By looking at the ansatz (2.4), we see that only two components of the vector field are non 
zero, and A^; where the subscripts refer to spatial indices and superscripts are internal indices. 
Projecting this vector on the surface, by stereographic projection, we see that the profile of the 
vector is that of a vortex. These vectors do not have definite charge under cr''; but, after some 
simple algebra we can change the basis to have definite charge vectors. As usual we can define 

= lb ir"^ and write 

^ = — dt + rr+ e+ + r" e_l (3.11) 

er e 

where the e± are: 

e± = sin{e) dip ± i dO. (3.12) 

As it is clear, by multiplying these basis with a length factor R, the real part of the basis is 
the azimuthal and the imaginary part of the basis is radial differential length elements in the 
complex plane. In figure (1) we have shown the directed curves for basis. At the poles the 
basis is completely polar, and at the equator there are equal contributions to the e±, from 9 and ip 
coordinates. The profile of the vector in this basis is also similar to a vortex. In hedgehog gauge 
the magnetic charge is related to 7r2(5'^), but, in abelian gauge it is the magnetic flow which is 
carried by the Dirac string [l2]; since the Dirac string passes through the center of the vortex, so 
the magnetic charge is the magnetic flow which is trapped by the vortex. In figure (2), we have 
shown the profile of the vector condensate, in hedgehog gauge. 

4 Numerical calculations 

In this section we construct a numerical self-garvitating Julia-Zee dyon, and study the behavior of 
the v.e.vs as functions of the sources of the model. The equation of motion for the matter functions 
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are: 



Figure 2: profile of tlie vector condensate. 



j"-^(x + yy(^)'-e-^ = o, 

K" + -{X- Y)'K' - e^'^iK^ + - e"^ - l) 
2 r 



0, 



(4.1) 



and tlie equations for metric functions are: 



-Y 



'"^-(ry'-i) + i,+ ^ 



SvrG, 



L2 



=2 



-^-(rA + 1) + — + — = — 



L2 



=2 



[X" + \{X'f - \X'Y' + \{X - Y)'] + A = ^T,^ 



SttG 



^2 V 



TV' 



(4.2) 



where 



/^2 , V ^ _ [(_yj2^ 



.4 +V[(7)']' + 



+ eV(</.) 



-r 



(K2_i)2 e-(^+^)r,J,, 



2r4 



+ 



e-^^ ^[(-)] + 



2 r 



(4.3) 
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To find a solution we assume tliat there is a horizon at ru, and consider a solution with the near 
horizon expansion of the form: 

J(r) = Ji(r - ru) + J2{r - ruf + ■■■ 
H{r) = Ho + Hi{r - m) + H2{r - rnf + ... 

K{r) = Ko + Ki{r- ru) + i^2(r - ruf + ... (4.4) 
e^W =Xi{r- m) + X2{r- ruf + ... 
g-nr) =Yi{r- m) + Y2{r - mf + ... 

using the equations of motion it is easy to see that only Ji, Hq, Kq, Xi are independent coefficients, 
and other coefficients can be determined from them. 
Our system of equations has a scaling symmetry: 

(r, rn, t,L) ^ a (r, rn, t, L) 

o 1 2 (4-5) 

{e,ip,J,K,H,e^,e^,X,e) {e,ip, J, K, H,e^ , X, e) 

Thus, choosing a = 1/r//, we can set r// = 1; we then defin^ 

r 

X = — 

L = — (4.6) 

_ t 

(L in this definition should not be confused with L in (2.7)). So that by starting from these 
expansions and running to large radii, from our numerical calculations, we see that the asymptotic 
form of the functions are: 



J{x) 


= px - p + 




H{x) 


= Ho + ^ 

x 


+ ... 


K{x) 


x 


+ ... 




= = I3 + - 






= j5 + ... 






= -ly- ds^ ~ 


L2 



(4.7) 



dr"^ + ^ dx"^ + x^ dnl. 

Note that generally in our numerical calculations, c ^ 1. In order for the asymptotic speed of light 
to be 1, we have to rescale the time. In fact there is another scale symmetry in the equations of 
motion: 

X{r) ^ X{r) + 2/3 ; t ^ e'^i ; At ef^Af (4.8) 



*In fact by this rescaling all the quantities will be in the units of rn- 
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Using this scaling we can write 

r^t = CT. (4.9) 

So the correct physical time is 

f = th i = c th T. (4-10) 

Since what we calculate numerically are Qxt^Ott and A^, for calculating temperature, chemical 
potential and charge density (and all quantities which are related to time coordinate in someway) 
correctly, this rescaling must be considered in the calculations. So for the temperature we have: 

T = ^{VF'drV^g^f = ^-1_( ^5^^/=^ )^^^. (4.11) 

In the series considered below, we use in our numerical calculations, 



(4.12) 



g-- = Yi{x-l) + ... 
9tt = Xi{x - 1) + ... 

the correct temperature is: 

T = -^ ^/X^l■ (4.13) 

47r c rn 

Also concerning the chemical potential we have, 

Ar dr = ( Ar-^)df = Af df; (4.14) 
crH 

so that Af is the physical field. 

In order to calculate the vacuum expectation values, we should study the fluctuations on the 
background. As can be seen from (3.7) to get the background v.e.vs it is enough to consider the 
fluctuations of the same form as those of the background. So by varying the matter function^ 

H H h(r) 
> — + e 



er er e 



(4.15) 

er er e 
K K Q{r) 



and solving the linearized equations for the fluctuations on the background of the calculated nu- 
merical solution, with the same initial conditions as the background functions, we see that their 
asymptotic forms are the same as those of the background: 

h{r) + \ + ... 

P^r) = Po + ^ + ... (4.16) 

r 



Q(r)=Qo + — + ...; 

r 

so ii,Ho,l — Kq, are sources and p,Hi,Ki, are condensates; also ho,Po,Qo, are probes and 
hi, Pi,Qi, are responses. 



^For simplicity we consider only s-wave fluctuations 
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When the radius of the horizon is fixed, although we have four sources {T, fj., S^, Sa) in the 
model, the phase space is a three dimensional hypersurface. It is similar to equation (2.11), where 
temperature and chemical potential are constrained and the phase space is one dimensional. In 
this case we expect the constraint to be: 

= T + ^ - + bcp \- Sa — • (4-17) 

s s s s 

In figure (3) we have plotted some curves related to some numerical solutions in the (T, fi) plane, 
for fixed values of 5*0 and Sa- 

From the numerical calculations the following emerges: 

1. The solutions with the desired boundary conditions exist only for a limited range of 5"^ and 
Sa and also temperature (chemical potential). 

2. As is shown in figures (4,5), the vacuum expectation value of the vector operator generally 
has different signs in different regions of the phase space, and there is a two dimensional 
surface on which the vacuum expectation value is zero. The sign changes smoothly. This 
phenomenon can be interpreted vortex reversing 



3. At fi = S(j, = 0, the sign of {V) changes at {SA)cr = 0, or alternatively, (i^o)cr = 1- Since the 
radial nonabelian magnetic field is proportional to oc {1 — K{r)'^)f, the boundary magnetic 
field is {Br)boundary OC (1 — -K^q)^- So at fi = S^ = 0, the vortex reverses when the boundary 
magnetic field is reversed. 

4. As is clear from figures (4,5), at nonzero chemical potential, even when the boundary magnetic 
field is zero, the vacuum expectation value of the vector operator can be nonzero; So the vortex 
can be a condensate. In order to cancel the vortex, we should impose a boundary magnetic 
field. Increasing fi, increases (Ko)cr and alternatively {Br)boundary which is needed to cancel 
the vortex;and increasing S^, decreases {Ko)cr- 

5. When S^j, is small enough, similar to the vector operator, the vacuum expectation value of the 
scalar operator has different signs in different regions of the phase space, and there is a two 



''Note that, as is clear from the equations of motion, the absolute sign of the matter functions (J,H,K) is not 
important, although sign changing is meaningful. In this paper we have considered positive solutions, (J,H,K)> 0. 
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(a) 



(b) 



Figure 4: plot of the phase subspace when S/p = 0.1. The red region is where the v.e.vs are negative, 
and in black region they are positive. Plot (a) is for {V), and plot (b) is for (O) 




(a) (b) 

Figure 5: plot of the phase subspace when S,^ = 1.5. Plot (a) is for {V), and plot (b) is for (O). 
The negative region for (O) is disappeared. 

dimensional surface on which the vacuum expectation value is zero, but this time it occurs 
at low temperatures. This phenomenon can be related to an instability, because in [6], it 
is shown that such an event occurs for linear fluctuations which means that at some special 
value of parameters the equations admit a zero mode solution, as a marginal stable mode. 
For vector fluctuations this is not the case. 

5 Conclusion 

In this article we studied the vacuum expectation values of dual operators from Einstein- Yang- 
Mills-Higgs action in a holographic setup and found that the radial magnetic field which is the 
main property of the magnetic monopoles does not play explicit role in the vacuum expectation 
values. In particular we consider the Julia-Zee dyon, and saw that the v.e.v of dual scalar and 
vector operators can be nonzero, and the profile of the vector condensate is similar to a vortex. 
We also constructed a numerical self-gravitating dyon and studied the changing of the v.e.vs as a 
function of the sources. We observed a sign changing phenomenon in the v.e.vs by changing tem- 
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perature, chemical potential, and scalar and vector sources. It will be interesting to check whether 
these sign changings relate to an instability. Also it seems reasonable to expect that the multi 
monopole solutions which are studied in \9j, correspond to multi vortex configurations in the dual 
field theory side. Any firm statement on this requires knowledge of the profile of the vector field in 
the bulk. 
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